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S <r- Default sample size 



Optionally receive S <- Desired sample size 



Generate (Yj,Ij) <- (0,0), j=l,2,...,S 



Initialize a reservoir R to empty 



M 0 



Provide seed number for random number generator 
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U <r- Random number, 0 < U < 1 , 
uniformly distributed, non-repeating 



Find (Yk,Ik) where < Yj, 1 < j < S 
(Minimum Y) 




Skip next record 



M<-M+l 



Transfer next record to reservoir as Rm 



(Yk,Ik) ^ (U,M) 



Arrange (Y,I) in a heap with respect to Y 
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Fig. 2 



